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Abstract 

We perform a Bayesian analysis in the basis of a recently developed two-paranreter 
class of hybrid equations of state that allow for high-mass compact star twins. While 
recently a wide range of radii, from 9-15 km, has been inferred for different neutron 
stars using different techniques, we perform our analysis under the supposition that 
the radii are towards the large end (13—15 km). We use this radius constraint together 
with the undebated statistically independent constraint for high masses (~ 2 M 0 ) as 
priors in selecting the most probable hybrid equations of state from a family with two 
free parameters: the baryon excluded volume in the hadronic phase and the 8-quark 
vector channel interaction in the quark matter phase. 


1 Introduction 

The study of the internal composition of neutron stars is an active held of research 
which relies on astrophysical observations that allow to refine theoretical models. In 
this respect, the recent observations of massive neutron stars nua have imposed 
important constraints on the stiffness of the equation of state (EoS), and therefore 
on the density ranges covered by the density profiles of such high-mass compact star 
interiors. On the other hand radius measurements are still far from being precise, 
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with most of the methods relying on either indirect measurements or model dependent 
assumptions like, e.g., for the neutron star atmospheres. While there is a wide range 
of claimed radii starting from, e.g., ~ 9 km [3] to ~ 15 km [3] we will proceed here 
on the assumption that the actual radii are large, as reported by [SI E]- For a recent 
review of astrophysical constraints on dense matter see, e.g., Ref. [7], 

In this contribution, we present a Bayesian analysis (BA) study case with a class 
of hybrid EoS characterized by two parameters. The first one stands for the baryonic 
exluded volume that determines the stiffness of hadronic matter at high densities. The 
second one is the coupling strength for an 8-quark vector current interaction which 
regulates the stiffness of the high-density quark matter phase. It turns out that 
within the range of variation for the excluded volume parameter there is a qualitative 
change in the mass-radius relation for high-mass stars: beyond a certain value for 
the excluded volume the first oder phase transition proceeds with a sufficiently large 
jump in the energy density to cause an instability which, thanks to the stiffness of 
the quark matter phase at high densities goes over to a stable sequence of hybrid 
stars, the so-called ’’third family” of compact stars. In this situation, the conditions 
are fulfilled for the high-mass twin phenomenon, where stars on the high-mass end of 
the second family of purely hadronic neutron stars are degenerate in mass with their 
twin stars on the lower mass part of the third family of hybrid stars bearing a quark 
matter core, see [8j and references therein. For a recent classification of stable hybrid 
star sequences under generic conditions for the EoS, see Ref. [9j. We note that the 
high-mass twin phenomenon is quite substantially based on a stiffening of both, the 
hadronic and the quark matter EoS towards higher densities; it is not obtained within 
a systematic parameter scan of hadronic vs. NJL quark matter EoS ra which lacks 
additional stiffening effects as those introduced in Ref. |8]. 

This possibility of high-mass twin stars is of great importance not only due to 
the possibility of identification of a critical endpoint in the QCD phase diagram mi 
IT2] but also because it provides a solution to several issues discussed in [[13]: the 
hyperon puzzle [IT], the masquerade problem [T3] and the reconfinement case [fB;, 
E]. Moreover, the transition between twin stars bears an energy reservoir [T5] that 
qualifies it as a possible engine for most energetic explosive astrophysical phenomena 
like gamma-ray bursts and fast radio bursts or play a role in contributing to the 
complex mechanism of core collapse supernova explosions. The present BA in the 
restricted two-dimensional parameter space of the new hybrid EoS performed with 
the modern high mass and large radius priors will give an answer to the question 
whether the high-mass twin star case is preferable over the connected hybrid star 
branch alternative. 
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2 Hybrid compact stars: hadronic and quark EoS 


For the neutron star description we utilize the density dependent relativistic mean- 
field EoS (DD2 [19, 20]) with an excluded volume correction that takes into account 
the quark substructure of nucleons (protons and neutrons) leading to Pauli block¬ 
ing effects. This correction is applied at suprasaturation densities and has the effect 
of stiffening the EoS without modifying any of the experimentally well constrained 
properties below and around saturation. In order to quantify the excluded volume 
effects, we introduce the closest packing parameter v = 100 x n^/fm -3 , where n v 
is the closest packing density. Increasing the excluded volume leads to a lowering of 
the closest packing parameter and therefore to a lowering of the critical pressure and 
energy density where the first order phase transition to the high density quark matter 
phase occurs, implemented in the form of a Maxwell construction. 

This high-density phase is described by a NJL EoS with multiquark interactions 
as introduced in EH- In this model the 774 parameter (related to the vector channels) 
serves to stiffen the EoS at higher densities, an important effect that helps to support 
hybrid stars as massive as the 2 M 0 stars recently detected. 

As for the observational constraints we retain the mass and radius measurements 
while ignoring the baryonic mass constraint of Podsiadlowski et al. [ 22 ] or Kitaura 
et al. [23j that were included in our earlier Bayesian analyses of Refs. [211125]. 

In the BA by Steiner et al. [26] the luminosity radius extracted for burst sources 
has been used to constrain a combined mass-radius relationship. This method is 
problematic, in particular because of the unknown stellar atmosphere composition, 
uncertainties in the distance to the source, the bias of the parabolic M-R constraint 
with the shape of stellar sequences in the M-R diagram for typical EoS and last but 
not least due to unknown details of the burst mechanism. 


3 Bayesian analysis technique 


We start by defining a vector of free parameters If = {z/, 774 }, which correspond to 
all the possible models with phase transition from nuclear to quark matter using the 
EoS described above. The way we sample these parameters is 


iti = 


1 V4{1) 


}■ 


(i) 


where i — 0 ... N — 1 with N — A r \ x N 2 such that i = N 2 x k + 1 and k = 0 ... Ni — 1, 
l — 0 ... N 2 — 1, with Ni and N 2 being the total number of parameters vn^ and 774 m, 
respectively. 

Once all the models are defined the neutron star configurations are obtained by 
solving the Tolman-Oppenheimer-Volkoff (TOV) equations [27. 28], These results al¬ 
low us to use different neutron star observations in order to determine the probability 
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that a given EoS model fulfils the observational constraints. We use here a mass 
constraint [2] and a radius constraint [5], 

Our goal is to find the set of most probable it i matching the above constraints 
using the BA technique. For initializing the BA we propose that a priori each vector 
of parameters it i has the same probability, P (it i) = 1/1V, for all i. 


3.1 Mass constraint 


We describe the error for the event Ea of a mass measurement of the high-mass pulsar 
PSR J0348+0432 | 2 ] with a normal distribution J\T(pa, <r\), where the mean value of 
the mass is pa = 2.01 M 0 and the variance is a a = 0.04 M 0 . Using this assumption 
we compute the conditional probability of the event Ea (under the condition that the 
neutron star is described by the EoS model with the parameters If i) with 

P(E a \^ i) = ^(Mi,p A ,a A ) ■ (2) 

Here Mj is the maximum mass accessible with the vector it t and &(x,p,a) is the 
cumulative distribution function for the Gaussian distribution 




1 

2 


1 + erf 


x — p\ 

V2^J 


( 3 ) 


3.2 Radius constraint 

We consider here a very promising technique to measure radii of neutron stars that 
is based on the pulse phase resolved X-ray spectroscopy which properly accounts for 
the system geometry of a radio pulsar. This radius measurement gives pb = 15.5 km 
and <jb = 1-5 km for PSR J0437-4715 [5]. Moreover, Hambaryan et al. [6] have also 
reported compatible radius measurements for RXJ 1856.5-3754. 

We compute the conditional probability of the event Eb that the measured radius 
of the neutron star corresponds to the model with If j as 

P (E b 1 7 ^ i) = $(Ri, Pb, &b) • ( 4 ) 

Here the value Ri is the maximum radius for the given vector If j. 


3.3 Calculation of a posteriori probabilities 

It is important to note that these measurements are independent of each other. This 
means that we can compute the complete conditional probability of an event E given 
iti that corresponds to the product of the conditional probabilities of all measure¬ 
ments, in our case resulting from the constraints Ea, Eb, 

P {E |*b) = P {E a ■ P {E b ). (5) 
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Thus, we can derive the probability of the measurement of an EoS represented by a 
vector of parameters it t using Bayes’ theorem 


3=0 


( 6 ) 


Since all P(lti ) = 1/N and P(lti\E) is normalized to P(E\ltj) = N, we have 

P{-ti\E) = P{E\-ti)/N . (7) 



Figure 1: Family of hybrid EoS (black lines) with a first order phase transition be¬ 
tween hadronic matter described by the excluded volume corrected DD2 relativistic 
meanfield EoS and quark matter described by a NJL model with higher order ( 8 - 
quark) interactions that provide the important stiffening at high densities due to a 
large coupling strength r /4 in the vector current channel, see [ 8 ] - The optimal EoS 
(bold red line) corresponds to an available volume parameter v = 33 and rj 4 = 5. The 
limits of the la range are indicated by bold green lines. 
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4 Results and discussion 


In order to cover a relevant set of possible hybrid compact star EoS in the pressure- 
energy density plane we have varied the closest packing parameter in the range v = 
33,, 100 corresponding to closest packing densities of n v = 0.33, 0.34,..., 1.00 
firm 3 , and the 8 -quark vector current coupling in the range rj 4 = 0,1,2,..., 30. 

The hybrid EoS resulting from Maxwell constructions between all combinations 
of hadronic and quark matter EoS in this two-dimensional parameter space is shown 
in Fig. [l] 



__ 1 __ 1 _ mm _ 1 __ 

10 12 14 16 
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Figure 2: Hybrid compact star sequences described by the family of hybrid EoS in 
Fig.0 used for Bayesian studies together with observational high precision mass (la, 
2a and 3a horizontal bands) [2] and radius measurements (la, 2a and 3a vertical 
bands) [5]. 

Fig.0 shows the resulting compact star sequences in the mass-radius plane which 
have a vertical branch corresponding to hadronic stars and an almost horizontal 
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branch for hybrid stars with a quark matter core. The parameters are chosen in such 
a range that their variation entails that the hybrid star branch becomes disconnected 
from the hadronic one due to the appearance of a set of unstable configurations. This 
characterizes the appearance of high-mass twin stars and happens in particular when 
varying the closest packing parameter v (i.e., the excluded volume) in all of these 
models. 

We have already observed in [ 8 ] that increasing the r / 4 parameter for a fixed v 
value corresponds to an increase of the maximum mass, while the difference between 
the radius of the hybrid star and its hadronic twin decreases. This general behaviour 
is preserved for any fixed v. 
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Figure 3: A posteriori probabilities derived from a two-dimensional Bayesian analysis. 
The vertical axis shows the probability for an EoS characterized by the parameters v 
and r /4 to fulfill the observational constraints for mass and radius. 

A posteriori probabilities of identification for these twin EoS are shown in Fig. [3j 
The models with higher probabilities correspond to the ones that predict large radii 


7 






and are associated to small closest packing parameters v. The optimal value of the 
vector coupling strength is 774 = 5. We note that the family of hybrid EoS accessible 
within the given two-dimensional parameter space describes sequences of compact 
stars that have a hadronic branch and a hybrid star branch. The latter appears 
connected to the former for not too large excluded volumina, i.e. for larger closest 
packing parameter v. For sufficiently small closest packing parameters the compact 
stars on the hadronic branch have radii exceeding 14 km and the phase transition 
proceeds with a sufficiently large jump in the energy density so that the hybrid star 
sequence gets disconnected from the hadronic one, forming a so-called ”third family”. 
We summarize our conclusions from this study: 

1. The most probable models exhibit high-mass twin star configurations with quite 
distinguishable radii, differing by about 2 km. 

2. The region of the most probable models in the two-dimensional parameter space 
is sufficiently narrow, covering the ranges 33 < v < 38 and 3 < 774 < 7. 

3. The most probable models have a relatively small closest packing parameter 
v ~ 33 and a not too large vector coupling strength 774 ~ 5. 

4. The existence of the horizontal branch signals a strong first order deconfinement 
phase transition and is a feature accessible to verification by observation. To 
that end, at least for two high-mass pulsars with masses ~ 2 M 0 (like PSR 
J1614-2230 and PSR J0348+0432) the radii should be measured to sufficient 
accuracy and turn out to be significantly different. 

The next two steps in the development of the approach are devoted to an improve¬ 
ment of the variability of the dense matter EoS within a two-dimensional parameter 
space embodying, e.g., also the purely hadronic case without a phase transition and 
to mimicking the occurrence of structures (so-called ’’pasta phases”) in the phase 
transition region [29j [30]. 
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